Nonequilibrium charge transport in quantum SINIS structures 
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Charge transport in a high-transmission single-mode long SINIS junction (S stands for supercon- 
ductor, I is an insulator, and N is a normal metal) is considered in the limit of low bias voltages and 
low temperatures. The kinetic equation for the quasiparticle distribution on the Andreev levels is 
derived taking into account both inelastic relaxation and volt age- driven Zener transitions between 
the levels. We show that, for a long junction when the number of levels is large, the Zener transitions 
enhance the action of each other and lead to a drastic increase of the dc current far above the critical 
Josephson current of the junction. 

PACS numbers: 73.23.-b, 74.45. +c, 74.78.Na 



I. INTRODUCTION 

Weak links consisting of superconducting (S) and nor- 
mal metal (N) parts separated by insulators (I) are the 
subject of intensive experimental and theoretical studies. 
Transport through these systems is determined by several 
factors, such as transparency of the contacts, specifics 
of the weak link area and energy relaxation in the junc- 
tion, leading to a series of nonlinear characteristics in the 
current-voltage relation. In stationary regime, numerous 
configurations have been analyzed (see Ref. [l[ for re- 
view), revealing the importance of both Andree\*2 and 
normal reflection processes taking place in the contact. 

When the injection rate of new particles into the junc- 
tion is greater than the corresponding rate of inelastic 
relaxation the nonequilibrium effects are to be taken into 
account. In realistic junctions with high current den- 
sity and especially at low temperatures the inelastic re- 
laxation may become less effective, resulting in a strong 
nonequilibrium which crucially affects the current trans- 
ported through the contact. For even a small bias voltage 
below the energy gap, the oscillating Josephson supercur- 
rent may be accompanied with a nonzero dc component 
corresponding to the dissipative processes. Nonequilib- 
rium situations in various SINIS-type junctions ranging 
from a point contact to a ballistic junction with finite 
length have been studied by many authors, the discus- 
sion to some extent has also concerned 3 -^ the inelastic 
relaxation effects. As is well known, the dc component 
exhibits a subgap structure at bias voltages eV = |A|/n 
which is associated with multiple Andreev reflections 
(MAR) 3 - 6,7 . The quasiparticles trapped in the junction 
are accelerated by the applied voltage, while, for each 
cycle of repeated electron-hole reflections at the two NS 
interfaces, the energy of the particle increases by 2eV 
until the accumulated energy enables it to escape the 
pair-potential well. This works in a broad voltage range, 
but becomes more and more complicated when relax- 
ation effects are included or the transparency at the in- 
terfaces differs from unity. The low-voltage MAR pro- 
cess for eV <C A in a ballistic contact is equivalent 
to the spectral flow along the Andreev energy levels^ 



where the phase difference (j) adiabatically depends on 
time <f> = 2eV/h+ <pQ. However, for a non-ideal trans- 
parency, the energy levels are separated from each other 
by minigaps (see the next section) which suppress the 
transitions from one level to the next thus cutting the 
spectral flow off. As a result, for very low voltages, the 
dc current is small for contacts with any realistic trans- 
parency T =^ 1. 

The interlevel transitions can take place by means of 
Zener tunneling near the avoided crossings of the An- 
dreev levels; they restore the spectral flow and give rise 
to a finite dissipative dc current. The Zener processes 
are more simple in short junctions (point contacts) where 
only two Andreev states corresponding to particles trav- 
elling in opposite directions exist; these levels have only 
one minigap at the phase difference <f) — ir. Effects of 
Zener tunneling on the transport properties of quantum 
point contacts have been studied in Refs. [HH; the dc 
current was found to have an exponential dependence on 
voltage in the low- voltage limit. 

For junctions where the center island has a length d 
longer than the superconducting coherence length £ the 
number of levels is proportional to the ratio g?/£ and can, 
thus, become large. If the transparency is not exactly 
unity, these levels are separated by minigaps at ip = nk, 
where k is an integer. In practice, such SINIS struc- 
tures can be made of a carbon nanotube or semiconduc- 
tor nanowire placed between two superconductors as in 
Refs. [9||lfj|. In Ref. [llj this type of junctions was sug- 
gested as a realization of a quantum charge pump where 
the minigaps were manipulated by the gate voltage be- 
ing sequentially closed in resonance with the Josephson 
frequency. In the present paper we consider the low- 
temperature charge transport in these junctions for con- 
stant bias and gate voltages. We derive the effective ki- 
netic equation for the quasiparticle distribution on the 
Andreev levels taking into account both the inelastic re- 
laxation on each level and the Zener transitions between 
the neighboring levels and demonstrate that the voltage- 
driven Zener transitions from one level to the next en- 
hance the action of each other and lead to a drastic 
increase of the dc current as the transition probability 
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grows with the applied voltage. 

We begin with a brief description of the spectral prop- 
erties of double-barrier SINIS structures in Section [ill 
In Section IIIII we derive the kinetic equation that deter- 
mines the distribution function on the Andreev levels in 
the presence of Zener transitions and inelastic relaxation. 
In Sections llVI and IVl we calculate the dc current and dis- 
cuss the results. 



II. MODEL 

We consider a quantum SINIS contact consisting of two 
superconducting leads connected by a normal conductor 
that has a single conducting mode. The insulating bar- 
riers have a high transparency such that the contact is 
nearly ballistic. In this Section we briefly summarize the 
spectral properties of SINIS contacts which are important 
for the transport characteristics. It is well known that the 
supercurrent flowing through such contact is determined 
by the Andreev states formed in the normal conductor 
and extended into the superconducting leads. The states 
can be described by the Bogoliubov-deGennes equations 
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E F + U{x) 



<5>V> + Hij) = eip , (1) 



where a z is the Pauli matrix in Nambu space, and 
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The superconducting gap is A = |A|e ±!(; °' /2 for x > d/2 
and x < —d/2, respectively, while A = for —d/2 < x < 
d/2. For simplicity we model the normal reflections at 
the interfaces as being produced by 5-function barriers 
U(x) = I6(x - d/2) + I5(x + d/2) assuming that the 
quasiparticle velocity in the superconducting leads is the 
same as in the normal conductor. 

In the normal region the particle, e ±lq+x , and hole, 
e ^iq-x ^ waves have amplitudes u and v^, respectively. 
The upper or lower signs refer to the waves propagating 
to the right = (u + , v~) or to the left ip < — (u~ , v + ). 
The particle (hole) momentum is q± — k x ± e/hv x where 
v x is the quasiparticle velocity of the mode and k x = 
mv x /h. Scattering at the right and left barriers couples 
the amplitudes of incident and reflected wavesi^: 



S R $> 



R- 
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(2) 

The scattering matrices for the right and left barriers 
are S R = S(xr) an d S L — S(xl), respectively, where 
Xl = ^0/2 while xr = 4>l^- The scattering matrices are 
unitary S'S = 1 because of conservation of the quasipar- 
ticle flux. Components of the S matrix for <5-like barriers 
and energies |e| < |A| are£ 



(t/ 2 -y 2 )|z|V^TT 

N U 2 + (U 2 - V 2 )Z 2 ' A 



uv 



IJ2 + (ip _ v 2 )Z 2 



Here Z — ml/h 2 k x 
2-V2[i + iA 



is the barrier strength and U 



V = U*. The scattering 



phase S is introduced through cotd = Z. Applying the 
scattering conditions at both ends of the normal region 
one can derive a compact equation for the spectrum of a 
SINIS contacOi 



ISVpsin 2 a' + |SU| 2 cos 2 (^/2) 



sin 2 (/3 + 7) 



(3) 



Here a = k x d, a' = a + S, and [3 = ed/hv x . The phase 
7 is defined as Sn = e ,,y |iSjv|. For short contacts d <g; 
hv x /\A\, the spectrum has two branches varying from 
e = ±|A|at<jf> = and separated from each other by 
minigaps at = it + 2irk. The energy spectra of SINIS 
contacts in various limits have been extensively studied 
by many authors^. 
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FIG. 1: (Color online) Examples of the spectra, Eq. ([3]), for 
a long SINIS contact with Z = 0.5 and \A\d/hv = 10. Dark 
(black online) lines: resonance sin a' = 0, the gaps disap- 
pear for <f) = 7r + 2nk; light (red online) lines: anti-resonance 
cos a' = 0, the gaps disappear for <f> = 2nk. 

In what follows we focus on long contacts, d ^> hv x /\A\ 
which have a large number of levels N ~ d\A\/hv x . These 
levels split off from the states with e = ±|A| and fill the 
energy interval — |A| < e < |A| with spacings of the 
order of Hv x /d. Examples of the spectra are shown in 
Fig. [TJ Each level is a function of the phase difference 
(j>; its range of variation is of the order of the interlevel 
spacing. The levels approach each other more closely 
at 4> = 7rfc where they are separated by minigaps. All 
minigaps at <f> — n(l + 2k) disappear for the resonance 
condition, sin a' = 0. Similarly, all minigaps at <j> = 2nk 
disappear for anti-resonance, |sina'| = 1. This follows 
from Eq. © due to unitarity ISWI 2 + IS^ 2 = 1. The 
low-energy levels, ei |A|, 7 -c 1 have the form 



ei = ±e + nhv x l/d 



(4) 



where 



hvr, 



eo = —j- arcsin y T 2 cos 2 — + (1 — T 2 ) sin 2 a' , (5) 

I is an integer, T = (1 + 2Z 2 )~ 1 is the transmission coef- 
ficient of the contact. The energy gaps at <j> = tt(1 + 2k) 
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are all equal, 



Se v = (2hv x /d) arcsin | sin a' \ \/\ — T 2 



The gaps btiix at = 2nk are given by the same ex- 
pression where | sin a' | is replaced with | cos a' \ . For a 
transparent contact, T = 1, all minigaps disappear. 



III. KINETIC EQUATION 

If a bias voltage V is applied across the superconduct- 
ing leads, the current through the contact has both ac 
and dc components. For low voltages, eV much smaller 
than A, the dc current is small for contacts with any 
transparency T 1. This is due to the presence of mini- 
gaps discussed in the previous Section. In long contacts 
the minigaps exist at (p = -nk and suppress the transitions 
between the levels thus preserving the equilibrium distri- 
bution of excitations. As a result, the current through 
the contact is simply the equilibrium supercurrent with 
the phase difference <p adiabatically depending on time, 
cj) = ujjt + <po where ujj = 2eV/H is the Josephson fre- 
quency. The dc component should thus vanish. However, 
the time dependence of the phase induces the Zener tran- 
sitions between the levels near the avoided crossing points 
at <j) = 7rfc which produce deviation from equilibrium. 
Since the transparency of the contact is of the order of 
unity, T ~ 1, particles with |e| > | A| have enough time to 
escape from the double-barrier region and to relax in the 
continuum. Particles at the continuum edges with en- 
ergies e = ±|A|, which are in equilibrium with the heat 
bath, are captured on the outermost Andreev levels when 
the latter split off from the continuum as the phase varies 
in time. Next these particles are excited to the neighbor- 
ing levels due to the interlevel transitions. Relaxation of 
thus created nonequilibrium distribution gives rise to a 
finite dc current. 

In this Section we derive the effective kinetic equation 
that describes the distribution function on the Andreev 
levels taking into account both the inelastic relaxation on 
each level and the Zener transitions between the neigh- 
boring levels. We assume that the Zener transitions take 
place only near the avoided crossings at <p = ixk. This 
can be realized if two conditions are fulfilled. First, the 
minigaps Se^ and 5e2K should be much smaller than the 
average distance hv x /d between the levels, implying that 
the contact is almost ballistic with a transparency close 
to unity, 1 — T <C 1. Second, the applied voltage should 
be small eV <C hv x /d such that it cannot excite transi- 
tions between levels far from the avoided crossings. We 
also assume that, between the Zener tunneling events, 
the distribution function relaxes according to 

dfn _ fn fn 

~dt ~ t 
where r is an inelastic relaxation time which we assume 



equilibrium distribution. Here n labels the levels con- 
secutively from the lowermost level at e = — |A| up to 
the uppermost level at e = +|A|. We will consider low 
temperatures T <C Kv x /d such that 



/(°)-sign(6„) 



(6) 



Since the superconducting phase difference depends lin- 
early on time the distribution function can be written as 

fn{<t>) = fn ] + fn where 



= f n e- h ^ 2eVr 



(7) 



The amplitudes /„ of the decaying parts are to be found 
using the conditions at the transition points. 
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FIG. 2: Scheme of the energy levels as functions of <f>. The 
arrows show the direction of the spectral flow. Enlarged are 
shown avoided crossings of levels at (f> — 2nk and (f> — n + 2-Kk. 

We assume that the interlevel transitions at (f> = 2i:k 
have a (constant) tunneling probability po while the tran- 
sitions at (p = 7T + 2ir k have a probability p„ . We will see 
later that the assumption of constant probabilities is well 
justified. We denote the points 4> = 27r/c before and after 
the tunneling events as 0=F, respectively, while the points 
4> = it + 2ixk before and after the tunneling events are 
denoted as 7rq=. All the respective points 0=F on a level 
n are equivalent due to the 2n periodicity, so are all the 
points 7r=p. Using the evolution equation ([7]) we couple 
the distribution functions at the consecutive instants of 
the tunneling events (see Fig. ^ 



constant, and /„ 



(0) 



1 



2n 



(0) 



tanh(e„/2T) is the 



where v = nh/2eVT. Here we introduce the upper + (or 
— ) indices to indicate explicitly the distributions at the 
spectrum branches increasing (or decreasing) as functions 
of <p. The coefficients 4>^ and \n are the- amplitudes 
/„ in Eq. ([7]) defined for the intervals < cj) < ir and 
7r < <j) < 2-7T, respectively. 
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The tunneling events impose the relations 

/+(0+) = po/+_ x (0-) + (1 -po)/n (0-) , (8) 
/-(tt+) = p 7r /- +1 ( 7r -) + (l-p 7r )/+( 7r -) , (9) 

/ I 7+i(o+) = po/,; +2 (o-) + (i-po)/+ +1 (o-) ,(io) 

fn+l(«+) = P*#(*-) + (l-jO/n +l (*-). (11) 

illustrated in Fig. [5] For the coefficients ^ and \ they 
become 



where 



+boX^-i + (l-Po)x^]e- I/ , (12) 
+ [PoX^+ 2 + (1 -Po)x^ + i]e _ ", (13) 



and 

X^ = ftr[/Si - / f l 0) ] + KVWi + (1 -^+]e-", (14) 

X++i = PnlfP - ifi] + K^t + (1 -P T )^i]e-".(lB) 

According to our picture of the spectral flow through 
the Andreev levels, the boundary conditions are imposed 
at the continuum edges in such a way that, for the levels 
increasing as functions of <p, the distribution / + coin- 
cides with the equilibrium at e = — |A|, while, for de- 
creasing levels, the distribution /~ coincides with the 
equilibrium at e = +|A|. Since the bias voltage is low, 
eV <C hv x /d, the equilibrium function in both supercon- 
ducting electrodes can be taken as /'°- ) = tanh(e n /2T). 
Let us assume that trapping of particles from the contin- 
uum occurs at <f> — 0; the boundary conditions are then 
formulated for the function at e = — |A| and for ip~ 
ate = +IAI: 



|A| 



, C= 



+|A! 



= 



(16) 



In this case it is convenient to exclude the functions x 
using Eqs. and HU and solve Eqs. ([12]), fTg]) for the 
functions ip. 

We choose the level index n in such a way that e„ > 
for n > 1 and e„ < for n < 0. Equations (JTSJ) 
and (fl~3l) then couple the levels n and n ± 2. Since the 
temperature is low and the distribution is given by Eq. 
©, the r.h.s. of these equations vanish for all n ^ 0,2. 
Therefore, the coefficients for n > 2 and n < satisfy 
the homogeneous equations. We assume that there are 
N + 1 levels with positive energies and N + 1 levels with 
negative energies such that the outermost levels touch the 
continuum. Therefore, for the uppermost level n = N + 
1, the solution of Eqs. (fl2|) . (fl~3| satisfies the condition 
ippf_L_i — 0. Similarly, for the lowermost level n = — N 
the solution satisfies V'ijy = ^' With these conditions, 
the solutions for n > 2 are 



e _ r(JV _„ +1)w _ _ e r{N-n+X) w+ 



c > e -r(N-n+l) _ e r(N-n+l) 



W± 



p e ±r + p v e Tr 



2p p^ cosh(r) 



C + PaPnO- - e ±2r ) + P0 + Pir - ^PoPtt 

The solutions for n < have the form 



ib + 

rn 



p r{N+n) e -r(N+n) 



e r(N+n) w+ _ e -r(N+n) w _ 



(18) 



The effective relaxation rate r > is found from the 
determinant condition w + w^ — 1 which gives 

4poPtt(C + l)sinh 2 (r) = (C + p + p n - 2p p^) 2 

- (po +Pir - ZpoPir) 2 ■ (19) 

where £ = e 2l/ — 1. 

For an ideally transparent contact po = p v = 1 we 
find r = v. For strong relaxation, v 1, we also have 
r ~ v, and the distribution relaxes quickly. The most 
interesting limit for a general case Po,Ptt 7^ 1 is when 
inelastic relaxation is weak, v <C 1. We find in this limit 



sinh 2 (r) = v(y + p +p n ~ 2p p„) / 'p p v 



(20) 



The relaxation rate r can be either large or small de- 



pending on the probabilities. The inverse rate 



de- 



scribes the broadening of distribution over the energy 
states and plays the role of an effective temperature 
T c ff = r~ 1 (de/dn) — irhv x /2rd. The effective temper- 
ature can be much higher than the interlevel spacing if 
r < 1. 

The coefficients c > and c< are coupled through the so- 
lutions of four non- homogeneous equations resulting from 
two equations (|12[) and (jT3j) taken for two values n = 2 
and n = 0. Inspecting equations for other n we see that 
only the two coefficients tpQ and ip^ cannot be described 
by the solutions Eqs. (1 1 T[) and (|18p of the homogeneous 
equations. We write 



*i>t = vr > + , = vr < + , (21) 

V,+ = ^ +> + S> , V^ + = ^ +< + *o • (22) 
Here 6q X and are four new unknown coefficients. The 

coefficients V'f > an£ l V'd" > are defined to satisfy the homo- 
geneous equations for e n > and are given by Eq. (JT7J) ; 
the coefficients ipi < and ip^ satisfy the homogeneous 
equations for e n < and are given by Eq. (|18p . Inserting 
Eqs. (f2~lj) , (|22p into the four equations obtained for n = 2 
and n = from Eqs. (|12p and (|13p we find all the four 
coefficients 5. The result is ip^ — % I J \ > ■ V'o" = ? / ; cT < while 



K < -K > = er(fi 0) 

^ < -i>t > = e V {fi 
These two equations yield c > 



(0) 



/o (0) ) 



C [e rN (1 + e r w+) - e~ rN (l + e~ r W - 



(23) 
(24) 
C where 

2e u . (25) 



We put — = 2 for low temperatures. Therefore, 
(17) 

the distribution possesses the symmetry ip_ n = —ip n+ i- 



f(0) 
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IV. CURRENT 

The contribution to the current due to the deviation 
from equilibrium is 



/n 



2e x - den 7 

' h 2^ f)rh U 



(26) 



e„>0 



where f n = f n — . The sum runs only over the local- 
ized Andreev states because the continuum states relax 
quickly so that their distribution is almost in equilib- 
rium. The equilibrium supercurrent has been calculated 
in Ref. (see also Ref. [l[ for a review); it is an oscil- 
lating function of the phase difference and thus has no 
contribution to the dc current. 

Denote de^/dep increasing (decreasing) parts of the 
spectrum e n (<f>) as a function of </>. We have for the cur- 
rent averaged in time 




di 



n+l 



.(27) 



f+ + f?fn f- 
Jn+1 T d(j) Jn 



d<j> Jn+1 + 90 Jn i 

/ J n=2l + l 

The sum over I runs from to L = iV/2 where 
AT = 2dA/*n-hv x 



is the total number of levels with e n > 0, i.e., for both 
signs in Eq. Q • 

In Eq. (|2"T|) we can use the ballistic spectrum with T — > 
1. In this limit | S'jy- 1 = 0, = 1, thus the spectrum in 
Eq. takes the form 

cos(c/>/2) = ±sin(/3 + 7 ) . 

Calculating the energy derivative of this equation for long 
junctions d|A| 3> hv x , we find 

dcf> 2d ' 

We neglected dj/de compared to d/3/de which holds for 
all energy levels excluding those in a narrow region near 
the gap edge, 1 — |e|/|A| <§; (Hv x /d\A\) 2 . This means in 
fact that, neglecting this narrow region, we can use Eqs. 
(gD and © with T = 1 for all n. 
We have for < 6 < n 



fn = tfe 



± -h<t>/2eVT 



while for 7r < d> < 2ir 



± h(<p-x)/2eVT 

An c 



We obtain for v -C 1 



where i? 1 = e 2 /nh is the quantum of conductance, and 



*(P0,Ptt) 



JV/2 

E 



[(Vvt+i -i>n) + (xt - X n +l)]„= 



2/+1 



The combination of coefficients ^ and \ that enters 
the expression for the supercurrent for / > 1 can be writ- 
ten through the solutions Eq. (TP?) of the homogeneous 
equations (TT2"]) - (rT5|) . The term I = contains xt which 
is expressed through ipQ. However, one can check that 
the jump in ip^ from n = to n = 2 is compensated by 
the jump f[ 0) ~ in Eq. (fT5|), 

Consider the limit of low relaxation r < 1 provided 
Nr 3> 1. The limit r <C 1 is realized when ^ <C p^^p^. 
In this case Eqs. ([15]). (fT7| . and (J2SJ) give 



We have from Eq. ([29]) 



4i/ 



(29) 



Av N/2 

®(Po,Pt<) = — E' 



-r(2fc+l) 



2poPti 



(30) 

We keep ^ in the denominator since the combination pg + 
Pk — 2pop v vanishes when po , p v — > 1 . 

When the inelastic relaxation rate is so small that 
Nv -C 1 the effective relaxation r can decrease such that 
Nr <C 1. Since the product /V(p + Pn — 2poPn) is gen- 
erally not small we find from Eqs. (I14p . (fT5| . (fT7| . and 
13 



- i> n ) + bit - Xn+l) 



N{po + Ptt - ZpoPtt) + 1 ' 



and 



${Po,P,r) = 



ZpoPir 



(PO +Ptt~ ZPoPtt) + 1/N 



(31) 



Eqs. dSO]) and (01]) go one into another for N ~ The 
exact expression for $ is found from Eqs. (TH| , (TT5]) , (1 1 7[) . 
and (|25[) . One can approximate the function <I>(po,?>7r) by 
an interpolation between Eqs. (|30| and (|31| in the form 



$(P0,Ptt) 



2poPt. 



(PO + Ptt - ZpoPn) + 2/3 



(32) 



where 2/3 w max (z/, /V -1 ) w f + /V -1 . 

The probability of Zener tunneling can be easily cal- 
culated for the spectrum in the form of Eqs. dm, dSJ) if 
(1 — T 2 ) <C 1. The phase difference is (f> — ujt + </>o- As 
a function of time, the distance between two neighboring 
levels for <b close to ir is 



Se = 

where t is small and 



Hv x Tu>j 



t 2 



r 2 =4sm 2 a'(l-T 2 )/T 2 iu 



J ■ 
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Probability of Zener tunneling is 



p n = exp 
where 



-— Im 

n 



Sedt 



= exp 



w . 2 / 

sin a 



uj Q = ttv x (1 - T 2 )/Td 



The distance between two levels for <f> close to 4> = and 
the corresponding probability of Zener tunneling po are 
given by the same expressions with sin a' replaced with 



cos a 



In Eq. (|32|) the term with (3 is only important when po 
and p n are close to unity. Therefore, one can write 



^{P0,Ptt) 



^exp 


u 


cosh 









cos(2a') 



1 + /3 



(33) 

When the bias voltage is low ujj/ujo < 1, such that 
z/ -C Po,Ptt -C 1, we have 



${Po,P7r) = 2exp ( [1 + |cos(2a') 



(34) 



The low-voltage part is exponential due to small Zener 
probabilities. The exponent exhibits strong oscillations 
as a function of a' which can be manipulated by varying 
the gate voltage. 

For higher bias voltages, u)j/u)q ^ 1, we have po,p v — > 
1 and 



$(po,Ptt) 



wo + wj/3 w + 7r/2r + wj/2A^ 



For a fully ballistic contact with po = Pn = F i-e., wo = 
our result agrees qualitatively with Ref. For these 
voltages we have two regimes. The FV curve is linear 
I = V/R as long as eV < A(7kj + ttH/t). The effective 
conductance is 



1 



1 



TTV x /d 



Rq luq + 7t/2t 



Inelastic relaxation can be neglected if ujq 7t/2t. In 
this case the conductance is much larger than the con- 
ductance quantum Rq/R = (1 — T 2 ) -1 . It is interesting 



to note that the effective conductance is independent of 
the gate voltage: it contains the sum of two functions 
in the exponents for po and p^, i.e., (uiq/uij) cos 2 a' and 
{(jJq/uij) sin 2 a' , which obviously is independent of a' . 

With increasing voltage up to eV > N(fku + nh/r) 
the I-V curve saturates at the value 

I = Nev x /d = 2e\A\/Trh 

which is by a factor N ^> 1 larger than the critical 
Josephson current of the junction 14 , I c ~ ev x /d. 

V. CONCLUSIONS 

To summarize, we have considered the low temperature 
charge transport in a nearly ballistic single-mode SINIS 
junction having a length d longer than the superconduct- 
ing coherence length £. In this junction the energy spec- 
trum of Andreev states has a large number of levels sep- 
arated from each other by minigaps which do not vanish 
in a realistic case when the transmission is not exactly 
unity. In the limit of low bias voltages, we have derived 
and solved the kinetic equation for the quasiparticle dis- 
tribution on the Andreev levels that takes into account 
both inelastic relaxation and voltage-driven Zener tran- 
sitions between the levels. We have shown that the Zener 
transitions enhance the action of each other and lead to 
a drastic increase of the dc current. The voltage depen- 
dence of the dc current is first exponential due to small 
probabilities of Zener tunneling. Next it goes over into a 
linear relation such that, at low temperatures when the 
inelastic relaxation rate is slow, its slope is determined 
by the average minigap in the spectrum. At higher volt- 
ages when the Zener probabilities approach unity, the 
dc current saturates at a value far exceeding the critical 
Josephson current of the junction. 
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